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НЕЕ. BALIE (Rhododendron spinuliferum) З: ЧЭ Ed P8 РА АХ БЕЗ НО SER Z5 НІЛ. ON ТЭ BR HE ТИЕ 
ЖЛЕ [Н] ЭЕ ВАЧ Е, WARAH FIASCO ОЗ АРЕВИЛЕЯ ТИЯСЫ ІШ Ж, Ж100 Xa A 
ПАН 28 Mit bid, HP 22 SAAS. ЖИБЕ ЛЕР ЛЕЛЕ 24 ARR 22 ^ GS ED ETT AY 
Wr, Жл: GAMARRA 2-5 EMEN, 3539289 3.4 4S, БОЯЙ A REA AR A BE 2190 7 0. 083 
~0. 792 ЖП 0. 153—0. 744, XTM AY 28 У DE de TE ARS AE EPI А ER Е Г К, AR ka 
A 22 ^C DB EID dE, ЖР 20 PASAT; ST ASI 2-6T x EDI, 3ESgg 3.27, Ж 
АДЛ ТАЯ BAe SPE 43-1 y 0. 000 ~ 0. 833 ЯП 0. 117-0. 736, JF E CT Ата n] ТЕЛЕН A 
MURIS a EEG ТЕ PT АП AR CE ЛЭХ 

жағы: MEERE; ЖЫЛЕ; ЕЖ; FIASCO; % УМ 
HATES: Q 943 ERIRE: A СЕ} =: 2095-0845 (2014) 01-041-06 
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Isolation and Characterization of Microsatellite Markers 
for the Chinese Endemic Species Rhododendron 
spinuliferum ( Ericaceae) 
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Kunming 650201, China; 3 University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: Rhododendron spinuliferum ( Ericaceae) is an ornamental and medicinal plant endemic to southwest Chi- 
na. In order to study hybridization between R. spinuliferum and R. spiciferum, a FIASCO (Fast Isolation by AFLP of 
Sequences Containing Repeats) method was used to develop microsatellite markers in R. spinuliferum. A total of 28 
microsatellite markers were isolated from 100 SSR primer pairs, of which 22 were polymorphic. Polymorphism of the 
22 polymorphic loci was assessed separately in 24 individuals collected from two wild populations. The number of al- 
leles per locus ranged from 2 to 5, with an average of 3. 4, while observed (Но) and expected (Н,) heterozygosi- 
ties varied from 0. 083 to 0. 792 and from 0. 153 to 0. 744, respectively. The same 28 microsatellite markers were al- 
so tested in two wild populations (12 individuals from each) of R. spiciferum. Twenty two of the markers were suc- 
cessfully amplified, of which 20 were polymorphic. Estimates of diversity in two natural populations of R. spiciferum 
based on the 20 polymorphic markers revealed that number of alleles per locus ranged from 2 to 6, with a mean of 


3.2, while observed (Ну) and expected ( H, ) heterozygosities ranged from 0. 000 to 0. 833 and from 0. 117 to 
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0. 736, respectively. These newly developed microsatellite markers will be used in future studies of hybridization and 


the population genetics of R. spinuliferum and its closely related species. 


Key words: Microsatellite markers; Rhododendron spinuliferum; Rhododendron spiciferum; FIASCO; Polymorphism 


Rhododendron L. is the largest genus in Ericace- 
ae including about 1 025 species, distributed from the 
northern temperate zones, throughout tropical south- 
eastern Asia to northeastern Australia ( Chamberlain 
et al., 1996). There are 571 species in China, of 
which 405 are endemic ( Fang et al., 2005). There is 
a long horticultural history of Rhododendron, and lots 
of species in this genus have been used as ornamental 
plants ( Yang et al., 1999). Rhododendron spinulife- 
rum Franch. is one of the important ornamental spe- 
cies, which is endemic to southwest China. The 
dried stems, leaves and flowers of R. spinuliferum 
can be used as Chinese folk medicine for eliminating 
phlegm, diminishing inflammation, relieving cough 
and asthma ( Chen et al., 1996). 

Natural hybridization has been recognized as an 
important factor of speciation and diversification 
within Ahododendron ( Milne et al., 2010). Numer- 
ous instances of natural hybridization have been re- 
ported in Rhododendron in previous studies ( Kron et 
al., 1993; Zhang et al., 2007; Milne and Abbott, 
2008; Ma et al., 2010; Zha et al., 2010). Rhodo- 
dendron X duclouxii is an inferred natural hybrid spe- 
cies between  R.spinuliferum апа R. spiciferum 
Franch. based on morphological characters ( Handel- 


Mazzetti, 1936). This had been confirmed by mo- 


lecular sequence data recently ( Yan et al., 2013). 


To further study the extent of hybridization between 
R. spinuliferum and R. spiciferum , codominant micro- 
satellite markers will be used. In this study, we tried 
to isolate and characterize suitable microsatellite 
markers from R. spinuliferum , and test their feasibili- 


ty in R. spiciferum samples. 


1 Materials and methods 
1.1 Materials 

Leaf samples of R. spinuliferum used to develop 
microsatellite markers were collected from two natu- 
ral populations in Yunnan province, China. Total of 
24 individuals from two R. spinuliferum populations 
(12 individuals from each) were selected to assess 
polymorphisms of the developed microsatellite mark- 
ers (Table 1). The feasibility of the developed mic- 
rosatellite markers was also assessed in 24 individu- 
als from two R. spiciferum natural populations ( 12 
individuals from each) (Table 1). Voucher speci- 
mens were deposited in the herbarium of the Kun- 
ming Institute of Botany, Chinese Academy of Sci- 
ences (KUN). 
1.2 Methods 
1.2.1 


Total genomic DNA was extracted from silica-gel- 


DNA extraction 


dried leaves using an improved CTAB ( cetyltrimethyl 
ammonium bromide) method (Liu and Gao, 2011). 


Table 1 Details of materials information used in this study 








Taxon Locality Gepgraphic Altitude/m Collection number 
SEEN NE EU YLJ-12008, YLJ-12009, YLJ-120010, YLJ-12011, 
R. spinuliferum B реш ы Ded m 2070 YLJ-12012, YLJ-12013, YLJ-12014, YLJ-12015, 
LA а : YLJ-12016, YLJ-12017, YLJ-12023, YLJ-12031 
Gen eae — YLJ-12802, YLJ-12803, YLJ-12804, YLJ-12805, 
К. spinuliferum LI : din е ЕТО? 1775 YLJ-12806, YLJ-12807, YLJ-12808, YLJ-12809, 
Бы фар ) YLJ-12810, YLJ-12811, YLJ-12813, YLJ-12817 
E E SR YIJ-12001, YLJ-12018, YLJ-12019, YLJ-12020, 
R. spiciferum Pune Mna л ман 2070 YLJ-12025, YLJ-12027, YLJ-12029, YLJ-12030, 
тена КЕ Нун ' YLJ-12036, YLJ-12037, YLJ-12043, YLJ-12051 
| YLJ-121089, YLJ-121093, YLJ-121097, YLJ-121111, 
- Yuxi, Eshan county, N 24.22? 
R. spiciferum à; 1830 Y1J-121120, YLJ-121122, YLJ-121123, YLJ-121124, 
Gaoping town E 102. 32? 


YLJ-121126, YLJ-121127, YLJ-121128, YLJ-121130 
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1.2.2  Microsatellite loci isolation, PCR amplifica- 
tion and polymorphism assessment 

The microsatellite loci were isolated based on 
the FIASCO protocol ( Zane et al., 2002). About 
300-500 ng genomic DNA was digested with Msel 
(New England Biolabs, Beverly, Massachusetts, USA) , 
and the digested DNA fragments were ligated to an 
Msel AFLP adaptor pair (5!-ТАСТСАССАСТСАТ- 
37/5 -САССАТСАСТССТСАС-3”) at 37 © for 2 h 
with T4 DNA ligase (Fermentas, Burlington, Ontar- 
io, Canada). Total of 5 uL of a diluted digestion- 
ligation mixture (1:10) was used for amplification 
reactions with the adaptor-specific primers Msel-N 
(5'-GATGAGTCCTGAGTAAN-3’) , with the follow- 
ing cycle program: 95 °C for 3 min, 30 cycles of 94 °C 
for 45 s, 50 °С for 60 s, 72 °C for 60 s, and a final 
extension step of 7 min at 72 ^C. The amplified frag- 
ments (200-800 bp) were enriched for microsatel- 
lite repeats by magnetic bead selection with 5'-bioti- 
nylated (АС), and (АС) ,;. These enriched frag- 
ments were amplified again with the MselI-N primers. 
The PCR products were purified using an EZNA Gel 
Extraction Kit ( Omega Bio-Tek, Guangzhou, Chi- 
na). The purified PCR products with enriched micro- 
satellite repeats were ligated into the pGEM-T vector 
(Promega, Madison, Wisconsin, USA) and trans- 
formed into DH5a cells (TaKaRa, Dalian, China). 
Identification of recombinant clones was performed in 
a blue/white selection assay. Positive clones were 
then tested for microsatellite inserts by PCR with 
(AC), (АС) and T7/Sp6 primers, respectively. 
Clones with positive inserts and appropriate size ( 300 
-700 bp) were then sequenced. The sequences of 
which contain microsatellite repeats (SSRs), and 
with sufficient flanking regions were then used for de- 
signing locus-specific primers with the program Oligo 
6. 0 (Offerman and Rychlik, 2003). 

The polymorphisms of all microsatellite loci 
were then assessed in 24 individuals from two natural 
populations (12 individuals from each) of R. spinu- 
liferum ( Table 1). PCR reactions were performed in 
20 wL volumes containing 50-100 ng genomic DNA, 


0. 6 uM of each primer, 7. 5 uL 2x Tag PCR Mas- 
terMix (containing 0. 1 U Taq polymerase/pL, 0.5 
mM dNTP each, 20 mM Tris-HCl (рН-8. 3), 100 
mM KCl, 3 mM MgCl, (Tiangen, Beijin, China). 
The PCR amplifications were conducted under the 
following conditions: 95 *C for 3 min followed by 32 
cycles at 94 % for 30s, an annealing temperature op- 
timized specifically for each primer pair ( Table 2) for 
45 s, 72°C for 60 s, and a final extension step at 72 °C 
for 7 min. The amplified fragments were separated on 
8% polyacrylamide denaturing gels with a 20 bp mo- 
lecular size standard ladder ( Fermentas, Burlington, 
Ontario, Canada) and visualized by silver staining. 
The polymorphic information content was calculated 
by PIC Calculator. Standard genetic diversity param- 
eters and deviations from the Hardy-Weinberg equi- 
librium were estimated in GENEPOP version 4. 0. 10 
( Rousset, 2008) for all polymorphic loci. Estimation 
for linkage disequilibrium between pairs of loci was 


performed also in GENEPOP version 4. 0. 10. 


2 Results and discussion 

A total of 294 clones with positive inserts and 
appropriate size were sequenced. Among these se- 
quences, 217 (73.896) sequences were found to 
contain microsatellite repeats ( SSRs) , and 133 of 
these sequences with sufficient flanking regions were 
suitable for designing locus-specific primers. Final- 
ly, total of 100 primer sets were designed for develo- 
ping microsatellite loci. 

The evaluation criteria for the amplification suc- 
cess rate of the loci followed Gao et al. (2012). Of 
the 100 primer pairs tested , 28 loci were successful- 
ly amplified, of which 22 showed polymorphisms, 
and six were monomorphic (Table 2). Sequences of 
these primers were deposited in GenBank under the 
accession numbers KC155596 to KC155623. For these 
polymorphic primers, the number of alleles per locus 
( A) ranged from 2 to 5, with an average of 3. 4, and 
the values for the observed (Но) and expected 
(Hj) heterozygosities ranged from 0. 083 to 0. 792 
and from 0. 153 to 0. 744, respectively. Five of the 
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22 polymorphic microsatellite loci deviated from Har- ed samples included. There was no significant geno- 
dy-Weinberg equilibrium (Р «0.01) (Table 3), typic linkage disequilibrium (LD) between any pair 
likely due to the presence of null alleles or few test- of loci at Р<0. 001. 


Table 2 Characteristics of 28 microsatellite loci developed in R. spinuliferum 


GenBank 


Locus Primer sequence (5'—3”) Repeat motif Size range/bp Та/°С . 
accession No. 





Polymorphic microsatellites 


F: TCTTCGTCTCCCTCTATCTTT 


Rh003 * ( TC)8 152-176 58 KC155596 
R; AACACACACAGACCTCAAATC 
Е; АТСАТТССТТСТТТТТСССТ 

Rh005 * (AG) 12 164-182 55 KC155597 
R; TCCACCCTCTGTCTCACTCT 
Е; TTGGAGTGAGAACAGAGAGG 

Rh008 * (AG) 14 202-234 55 KC155598 
R; TAATAGGCAGCATCTCCCAT 
F: GGTAGCCACACTGTTGAAAT 

Rh009 * (AG)8 216-230 54 KC155599 
R: CTTCCCCTCCATCTTGTTCT 
Е; TTTGGCTCATCGCTTTTAGT 

Rh017 * ( TC) 10 151-175 54 KC155600 
R: GAGAGCATCCAAGTCCCTAT 
F: GCATCTCAAGAACACAATA 

Rh020 * (AG)9 109-143 51 КС155601 
R; TCAAGAAGGTCCTCCCAGTC 
Ғ; GAGGAGAGAAAAGGACAAG 

Rh031 * (AG)14 231-237 49 KC155603 
К; AGTCTTCTTCCTTACCAACG 
F: GGGCAAACATTCATACATAA 

Rh032 * ( TC) 16 296-308 59 KC155604 
R; AGGCAGGCAGGCACCAGAAG 
Ғ; CAAAAAACACACCGCAGACG 


Rh034 * (AG)9 193-203 52 KC155605 
: ТСАТСССТССАТССАТААТ 


: CCTGGGCAAGAGAGAAAACT 
: ACAGCGATGGCGATTTGAAC 


: TCCTAATCCCTCCATCTCCC 
: GCCGTTCCATACAGTACCAA 


: CGATTGCCATTTGCCACTACCT 
Rh041 * ( TC)7 148-178 55 KC155608 
: CCACAACTCCGCTGCTACTG 


: AGTTCCCCAAATCTCTTCTC 

: TCATTTTCTTTTCTCTGCCT 

: TGTAGCAAACCCATCTCACC 

Rh054 (TC)8 261-275 58 KC155611 
: TCACCTGGGCATAACTAATC 


: GATATGGACTCCGACAAGGT 
Rh058 * ( TC)9 174-180 58 KC155612 


Rh037 (AG)8--- (AC) 11 279-287 55 KC155606 


Rh039 * (TC) 10 156-168 57 KC155607 


т mu a a a 


Rh043 * (А6)23 149-175 53 KC155610 


шн ш =ч m 


= ш 


R: GGCGAGATCGTGGAGAAAAT 
Ғ; AAGAGATTGGAAGGGTTGAT 

Rh060 * (АС)7--(АС)7 166-172 54 KC155613 
R; TCATAGTGTGGCAAAACGAC 
Е; TGACGACATGGGACTTTAGA 

Rh063 ( TC)20 164-172 32 KC155614 
R; ACCCTTTCTTCATCTTCCAG 
F: TAAAAAAATGGGGCTAAAGT 

Rh065 (AG) 16 261-283 50 KC155616 
R; GACATTGACGCAGCCGAACC 
Ғ; ССТСТАСССТТАТСАТТТТА 

Rh072* ( TC) 25 169-181 57 KC155617 
R; AAGACGGACGAAACACATC 
F: ATACACCACCATTCATACGC 

Rh076 * (AG)17G( AG)6---(AG)8 258-310 20 KC155618 
К; TAGAGAGTGGGGTTGATTAG 
Е; CAATGATGTGAAAGCCCTGG 

Rh078 * (TG)8 284-300 50 KC155619 
R; AGGATTCCAATTAGTAAACG 
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Table 2 continued 





GenBank 
Locus Primer sequence (5'-3') Repeat motif Size range/bp Та/°С | 
accession No. 

Е; ATCACCCAAGCAATAGTCTG 

Rh086 * (TC) 9-- (TG) 8 269-281 57 KC155620 
R; ATTTTCCACACGATACAGGC 

Monomorphic microsatellites 
Е; СТАССАТСААСАТСАСАСТС 

Rh023 * (TC)8 131 52 KC155602 
R; AGTAAAAAGAGAAGGGGAGT 
Е; САСААСТСТТССААСАТТСС 

Rh042 * ( TC) 7C( TC) 9 165 55 KC155609 
R; GACGGGAGTTATCGCTGAAG 
F: GATGGTAGTTTCAACGCAAG 

Rh064 * (AG)9 193 52 КС155615 
R; АСТССТТТСТТТТСТСАССТ 
Ғ; AGAATAGAAGGTTGAAGGGT 

Rh087 ( TC) 13 217 52 KC155621 
R; AAGGCTTGAATGAGGTTGAT 
Е: CCCTCCTCTCTCAACAAAAG 

Rh096 * (TC) 10 157 54 KC155622 
R; TCAGAGTTGTTCGGTGTGTG 
F: AAACCCCATTACAGTAGATT 

Rh098 (AG)9 189 50 KC155623 


К; ACTGGACCCTTGAAACCTAAC 


Note; Ta, PCR annealing temperature; " , successful amplification in R. spiciferum 


Table 3 Results of the polymorphic microsatellite loci evaluated in two wild populations ( 12 individuals from each) 


of R. spinuliferum and R. spicifeum respectively 


R. spinuliferum 


R. spicifeum 








Locus 

М, Ho Ay Puy Ny Ho Hy Puy 
Rh003 4. 000 0. 739 0. 616 0. 113 3. 000 0. 833 0. 542 0. 005 * 
Rh005 4. 000 0. 375 0. 318 1. 000 2. 000 0. 250 0.278 0. 502 
Rh008 4. 000 0. 417 0. 357 1. 000 4. 000 0. 667 0. 587 0. 142 
Rh009 2. 000 0. 083 0. 153 0. 128 2. 000 0. 545 0. 496 1. 000 
Rh017 5. 000 0. 750 0. 724 0. 254 2. 000 0. 500 0. 469 1. 000 
Rh020 4. 000 0. 792 0. 744 0. 090 2. 000 0. 542 0. 430 0. 355 
Rh023 2. 000 0. 042 0. 117 0. 063 
Rh031 4. 000 0. 333 0. 510 0. 011 4. 000 0. 500 0. 688 0. 148 
Rh032 3. 000 0. 565 0. 638 0. 336 4. 000 0. 333 0. 641 0. 001 * 
Rh034 3. 000 0. 250 0. 227 1. 000 4. 000 0. 708 0. 736 0. 023 
Rh037 4. 000 0. 227 0. 714 0. 000 * 
Rh039 4. 000 0. 667 0. 643 0. 565 3. 000 0. 542 0. 594 0. 477 
Rh041 4. 000 0. 583 0. 513 1. 000 3. 000 0. 500 0. 492 0.301 
Rh043 2. 000 0. 458 0. 430 1. 000 3. 000 0. 250 0. 624 0. 000 * 
Rh054 2. 000 0. 174 0. 476 0. 002 * 
Rh058 4. 000 0.417 0. 506 0. 091 4. 000 0. 792 0. 711 0. 104 
Rh060 3. 000 0. 250 0. 473 0. 001 * 3. 000 0. 333 0. 645 0. 008 * 
Rh063 3. 000 0. 625 0. 555 1. 000 
Rh065 3. 000 0. 708 0. 661 0. 670 
Rh072 3. 000 0. 208 0. 659 0. 000 * 5. 000 0. 708 0. 728 0. 523 
Rh076 4. 000 0. 375 0. 704 0. 000 * 6. 000 0. 333 0. 538 0. 016 
Rh078 3. 000 0. 333 0. 351 0. 157 3. 000 0. 333 0. 586 0. 015 
Rh086 3. 000 0. 208 0. 320 0. 056 2. 000 0. 000 0. 287 0. 000 * 
Rh096 2. 000 0. 208 0. 492 0. 006 * 


Note; №, , number of alleles revealed; Но, observed heterozygosity; Hj, expected heterozygosity ; 


from Hardy-Weinberg equilibrium ( P«0. 01). 


* 
, 


polymorphic microsatellite loci 


deviating 
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The 28 microsatellite markers were also tested 
іп R. spiciferum using the same PCR conditions as in 
R. spinuliferum. Of the 28 loci tested, 22 SSR mark- 
ers were amplified successfully, of which 20 loci 
showed polymorphisms and two loci were monomor- 
phic (Rh042 and Rh064) in R. spiciferum (Table 2 
& 3). The two monomorphic microsatellite markers 
( Rh023 and Rh096) in R. spinuliferum showed pol- 
ymorphism in R. spiciferum ( Table 2 & 3). For the 
20 polymorphic markers of R. spiciferum , the number 
of alleles per locus ranged from 2 to 6, with a mean 
of 3. 2. The observed (А) and expected ( H,) het- 
erozygosity ranged from 0.000 to 0.833 and from 
0. 117 to 0. 736, respectively. Six of the 20 polymor- 
phic microsatellite loci deviated from the Hardy- 
Weinberg equilibrium (P<0.01) (Table 3). 

In summary, of the 28 microsatellite markers 
firstly developed in R. spinuliferum, most worked in 
R. spiciferum (79%). Thus, these codominant mic- 
rosatellite markers developed in this study will be 
very useful to investigate the hybrid speciation sce- 
nario between R. spinuliferum and R. spiciferum , and 
also be useful to assess the genetic diversity and 
population structure of R. spinuliferum апа other 


closely related species. 
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